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Fluctuating hydrodynamics in a vertically vibrated granular fluid with gravity 

Giulio Costantini and Andrea Puglisi 
CNR-ISC and Dipartimento di Fisica, Umversitd Sapienza - p.le A. Mora 2, 00185, Roma, Italy 

We investigate hydrodynamic fluctuations in a 2D granular fluid excited by a vibrating base and 
in the presence of gravity, focusing on the transverse velocity modes. Since the system is inhomo- 
geneous, we measure fluctuations in horizontal layers whose width is smaller than the characteristic 
hydrodynamic lengths: they can be considered as almost-homogeneous subsystems. The large time 
decay of autocorrelations of modes is exponential and compatible with vorticity diffusion due to 
shear viscosity, as in equilibrium fluids. The velocity structure factor, which strongly deviates from 
the equilibrium constant behavior, is well reproduced by an effective fluctuating hydrodynamics 
described by two noise terms: the first associated with vorticity diffusion and the second with the 
local energy exchange, which have internal and external character, respectively. 

PACS numbers: 45.70.-n,05.40.-a 



I. INTRODUCTION 



(N 

Granular systems are ubiquitous on Earth: understanding those materials is essential for the optimization of a 
variety of industrial processes in the pharmaceutical, food, cosmetic, chemical, petroleum, polymer and ceramic 
^1^ industries [1, 2], it is also propaedeutic for modeling geophysical phenomena such as sediment fluidization in sudden 

£^ landslides [3-5] and soil mechanics [6]. Finally, it provides to be a challenge for the statistical physics of out-of- 
equilibrium complex materials [7-11]. Fluctuations in granular media are usually non-negligible because of their 
C^ ■ inherently "small" nature: they are usually composed of a few thousands particles, sometimes even less. On the 
other side equilibrium statistical mechanics, which is equipped with tools to describe large scale fluctuations of finite 
systems [12] , is not adequate for systems with non-conservative interactions, such as granular fluids undergoing inelastic 
'^ I collisions. A deep understanding of the single-particle velocity fluctuations has been obtained in recent years, through 
the study of the granular Boltzmann equation and the introduction of several effective models [13-16]. Fluctuations 
involving more than one degree of freedom, i.e. correlations, are instead less understood: the reason is that they 
not always have an equilibrium counterpart. For instance, velocity correlations are ruled out in homogeneous steady 
elastic fluids, while they play an important role in homogeneous steady inelastic gases [17-20]. Important progresses 
on the study of velocity correlations have been achieved for hydrodynamic (i.e. "slow") modes, in situations where a 
good separation of length and time scales is observed. While phenomenological fluctuating hydrodynamics has been 
^f-^ \ used in the first study of velocity correlations in homogeneous granular gases [17], a more rigorous treatment through 
\^ . projection formalism from microscopic models has been recently obtained [21]. One of the main results, seen both in 
^^ ' the theory and in simulations [22, 23], is that the hydrodynamic noise is not white and does not satisfy the second 
r — kind Fluctuation-Dissipation Relation (FDR). These two deviations from the equilibrium fluctuating hydrodynamics 
f^ ■ are observable but, however, relatively small in homogeneous systems. On the other side, a certain sensitivity to the 
particular modelling of energy injection has been observed [24, 25]. It is therefore tempting to extend this study to 
more realistic energy injection models, with the objective of providing a framework for future experiments and more 
formal analytical treatment. Here, in particular, we are interested in simulating a typical experimental setup, that 
is a 2D vertical box, under gravity, whose vibrating base acts as a thermal reservoir [24, 26-28]. The loss of kinetic 
energy, due to inelastic collisions in the bulk, is counterbalanced by the random floor vibration which injects energy 
to particles bouncing on it. This setup, which has received careful hydrodynamic treatment [26], is characterized by 
non-constant density and temperature profiles and by a steady (inhomogeneous) heat flux [29], carrying energy from 
the bottom region up to the top one. Our analysis, which is purely numerical, focuses on the velocity fields, which is 
zero on average, provided that convection rolls are not present in the system [30, 31]. Fluctuations of the shear modes 
(vertical velocity modulated in the horizontal direction at wavelength fc) are measured in subsystems small enough to 
be considered "homogeneous" : a subsystem here is a horizontal layer taking the full width of the system and a fraction 
of its height, smaller than the smaller characteristic length of hydrodynamic profiles. Three main characterizations 
of these fluctuations are in order: a) the large time decay of their autocorrelation {U±{k,t)U±{—k,0)) ^ exp[—q{k)t]; 
b) the dependence of {U±{k,0)U±{—k,0)) with fc, which is a constant for equilibrium fluids, and it is not in non- 
equilibrium ones; c) the relation between the two previous quantities, which at equilibrium is established in a simple 
form by the FDR. Previously, fluctuations have been studied in this particular "setup" focusing on other quantities. 
In [27] single particle velocity fluctuations and density correlations at different height have been studied. In [24] the 
fluctuations of the total kinetic energy have been investigated, showing that strong non-Gaussian behavior is not due 
to strong correlations but can be simply accounted by considering inhomogeneous temperature and local Gaussian 



velocity distributions. This consideration led us to the choice of studying sub-regions of the whole system where 
homogeneity can be assumed and velocity distributions are not far from the Gaussian. The drawback of this choice 
is, of course, that each sub-region is an open system and this fact must be taken into account. In Section II we define 
the model, the simulation scheme and the quantities under study. In Section III the separation of the inhomogeneous 
system into nearly homogeneous sub-systems is detailed. In Section IV we discuss the numerical results for the points 
a, b and c described above. In Section V we draw some conclusion and discuss perspectives. 

II. MODEL AND QUANTITIES UNDER STUDY 

We perform Molecular Dynamics simulations of inelastic hard disks of diameter a and mass 1, enclosed in a 2D 
rectangular volume of width Lx and height Ly, under the presence of gravity (pointing in the —y direction) with 
acceleration g. The energy is injected to the system by a thermal base: this means that a particle colliding with the 
wall at 2/ = is bounced back with its x velocity component Vx extracted by a Gaussian distribution with temperature 
Tft and with its y velocity component Vy extracted with the distribution [32] 
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P{vy)^^^e~^. (1) 

The ceiling of the box is elastically reflecting. Periodic boundary conditions in the x direction are adopted. The disks 
collide inelastically with usual rule dictating the instantaneous velocity change v^ — >■ v^ of particle i after a collision 
with particle j: 

1 + a 
Vj=v, —[{■Vi-v^)-n]h, (2) 

where a G [0, 1] is the restitution coefficient (= 1 for elastic systems) and h is the unit vector joining the colliding 
disks. We are interested in measuring the velocity shear mode, modulated along the x direction, which is the only 
one where periodic boundary conditions are enforced. This fluctuating variable with complex values is defined as 

N 

C/^(fc,t)=^«,,,(i)e-'=-^W. (3) 

where k is the wave number of chosen mode, N the total number of particles and Xj (t) is the x-coordinate of particle 

J at time t. 

The main function under investigation is the rescaled autocorrelation of mode k: 

c^^,^,^^iuAkS£l(Ml, (4) 

measured in the steady state, where w^^ — 2Tg^y and Tg^p — (v'i) is the "granular temperature" in the /3 direction, 
with /3 being x or y. Note that Tg^y/Tg^x ~ 1 in the region of interest (which includes most of the system, excluding 
the boundary layers near the bottom and top walls). We will also use the shorthand notation U±{t) = U±{kmm,t) 
and C±{t) = C±{kmin,t) for the largest mode kmin — ^tt/Lx- In particular we will inspect two main features of the 
above correlation function, i.e.: 



• 



• 



its large time decay, focusing on the dependence on the parameters such as the granular temperature Tg^y and 
the particle density n and on the wave number k; 

the variance of the fluctuations of the modes as a function of fc, C±{k, 0), which is nothing else than the rescaled 
velocity structure factor; that quantity gives information about the spatial extent of velocity correlations. 



Since the fluid receives energy from a boundary and gravity is present, all hydrodynamic profiles are inhomogeneous. 
In order to avoid the mixing of information coming from regions with different densities and temperatures, in the 
following section we discuss how to define sub-systems which can be assumed to be nearly homogeneous for our 
purposes. Before doing that, however, we discuss what is expected in previously studied homogeneous systems. 



1. Behavior of a homogeneous system 



For comparison, we briefly review the behavior of the above quantities in an equilibrium homogeneous 2D system 
and in a granular system both in the homogeneous cooling regime and in the homogeneous stationary state due to 
uniform random driving. 

At equilibrium (at temperature Tg — w^^/2), the Landau-Lifshitz fluctuating hydrodynamics, based on Einstein 
fluctuation formula, predicts a Langevin equation for the hydrodynamic shear modes of the kind 

dtUj_{k,t) = -vk^U^{k,t)+£,{k,t), 

where v is the kinematic viscosity and ^(t) is a white noise with zero average and 

(e(fe, mk'X)) = 5k',-k5{t - t')2TgNvk\ 

Based on this equation, one has 
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uk-'t 



(5) 
(6) 

(7) 



In the inelastic homogeneous cooling regime, which is also dilute, the amplitude of fluctuations is decaying: indeed, 
because of cooling, the kinematic viscosity v and Tg decrease with time. In these systems v can be written, in the 
first Sonine approximation, as 



v = g2{(T) 
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(8) 



where 52(c) is the equilibrium value of the pair-correlation function, A ~ [na)^^ is proportional to the mean free 
path and the coefficient 02(0) is a function of the restitution coefficient only (see [21] for details). 
For that reason one uses a new time-scale r which is proportional to the cumulative number of collisions [33] : under 
this new time-scale the system reaches a stationary regime. Rigorous treatment from a fluctuating Boltzmann equation 
in this regime yields [21] 



drU^{k,T) = - ( J.fc2 - % ) U^{k,T)+S,'{k,T) 



where 
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02(0) 



is the cooling rate (see [33] for a definition). The noise results to be non- white, with a correlation 

(e'(fc, T)i'{k', t')) - 4',-fe2iVTgfc2G(lr - r'l). 



(9) 



(10) 



(11) 



with G(lsl) 7^ 5{s) a function which is given in details in [21]. The steady autocorrelation C±{t) obtained from the 
above Langevin equation is not a simple exponential; anyway it has an exponential tail at large times. In particular 
one obtains 



Cx(fc,0) 



yj.fc2-Cff/2 

U2 



C.(fc,r).f%±^e-(--CW.). 
2 vk^ - QhI^ 
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(12) 
(13) 



where the two new coefficients vi and V2 (the latter is usually smaller than the first) are computed in [21]. In the 
elastic limit C,h ^ ^, vi ^ v and t^2 ~> and result (7) is recovered. In the inelastic case {C,h > 0) one immediately 
sees that the theory is limited to k large enough to have i^fc^ — Cff/2 > 0: modes with smaller wavenumbers are 
unstable. That condition obliges to consider systems smaller than a critical size to avoid the instability. 

Finally, we consider the case of a homogeneously driven granular gas, as described in [18], where the hydrodynamics 
of a gas of inelastic grains which receive energy by random uncorrelated velocity kicks is studied. In that work 
hydrodynamic fluctuations are described using an effective noise that is the sum of an internal and an external noise. 
The former is originated from the rapid fluctuations of microscopic degrees of freedom and its strength can be obtained 



from an FDR with respect to internal relaxation. The latter, instead, is due to the random accelerations received by 
particles from the external driving. The strength of this noise is such that, in the steady state, it balances the energy 
loss due to the collisions. The result for the shear mode with small inelasticity is a Langevin equation: 

drU±{k,t) ^-iyPU_L{k,t)+C'{k,t), (14) 

with 

(e"(fc, i)e"(fc', t')) = Sk',-kS(t - t')NTg {2ve + Ch) ■ (15) 

The corresponding autocorrelation C±{k,t) is written as 

We stress that other models exist to get a spatially homogeneous and stationary granular gas, for instance an important 
modification of the last one is that introduced in [16], where all particles feel the presence of a viscous bath: the 
fluctuating hydrodynamics for such a model has been recently studied [34] and demonstrated to fairly describe the 
experimental behavior of a quasi-2D system on a horizontal vibrating plate. 

An important property common to all the systems discussed above is the presence of spatial and thermal homogeneity 
that is, instead, absent in our system. We mention that fluctuating hydrodynamics in molecular (elastic) fluids where 
a non-equilibrium stationary state is obtained by imposing (temperature or density) gradients has been studied in 
the literature, see for instance [35, 36]. In particular, in [35] the authors analyze the fluctuations of a fluid under 
a stationary heat flux in the presence of a gravity fleld. They show that, for small temperature gradients, such 
system can be studied within a Landau-Lifshitz approach postulating that FDR still holds locally. This implies that 
the transverse velocity modes are well described by Eq. (5) where all parameters and coefficients take their local 
(position-dependent) value, and therefore do not display any correlation, yielding a flat structure factor. 

In the following we will see that the hydrodynamic fluctuations of the granular setup considered here, which is a 
nearly homogeneous sub-system belonging to a inhomogeneous one, are, for some aspects, well reproduced by a model 
similar to Eq. (14). 

III. QUASI-HOMOGENEOUS SUB-SYSTEMS 

The main external parameters influencing the regime are, besides the dimensions of the box, the number N of 
particles, the floor temperature Th and the restitution coefficient a. To give a flavour of the regime we are working 
on, in Fig. 1 we represent two typical configurations of particle positions in the system from simulations with different 
restitution coefficients. To fix ideas, useful for the following discussion, we now focus on a particular choice of 
parameters. The system has height Ly/a = 600 and width Lx/(J = 180, the number of particles is A^ = 1500 
(corresponding to a number of resting layers iV^ = Na/L^ « 8.3), and the restitution coefficient is a — 0.95. Gravity 
acceleration is set to g = 9.8, while Ti,/{gLy) — 2.55. 

In Figure 2 the behavior of the main "local" observables is reported versus the distance from the floor. We have 
divided the total height into mi = 10 horizontal layers of height Li/<t = 60, so that miLi — Ly. In each layer 
(whose height y is measured at its middle point) we have computed the average of some observable which is expected 
to be slowly varying, i.e. the number of particles in the layer N{y) and the average kinetic energy in that layer 
Tg{y) = ^[Tx{y) + Ty{y)], which are directly associated to the density and temperature hydrodynamic fields [37]. 
In Figures 2a and 2c we show N{y) and Tg{y) respectively. In Fig. 2a one can also appreciate (see right scale) 
the local packing fraction (j>{y) — n(y)7ro'^/4 where n(y) — N{y)/{LxLi) is the local density. We also give further 
information to assess the granular regime chosen for our study. In Figure 2b we report the behavior of the local 
mean free path X{y) = (n{y)a)~^ (note that X{y) is proportional to the local mean free path through a order 1 
geometrical factor and through the Enskog constant g2(c) which, in our dilute case, is very close to 1). Finally in 
Fig. 2d some local characteristic times are reported for the different layers: Tc{y) — \{y) / vth{y) is the mean collision 
time, Texitiy) = Li/vth{y) is the mean exit time from the layer with Vthiy) equal to the local thermal velocity and r^iy) 
is the time associated to vorticity diffusion obtained from the expectation for the local kinematic viscosity v{y), taken 
from the granular Enskog theory [37] . These times will be useful (and will be explained in detail) in the discussion of 
Section IV. 

We have verified that the choice of the system parameters as well as the choice of mi meet the following constraints: 

• Stability of horizontal modes: if the system is too large it is known that instabilities can arise, starting with shear 
modes and then involving clustering [17, 38]; a minimum criterion to avoid that is enforcing k^in > kj_{y) = 
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FIG. 1: Snapshot of the system with A'' — 1500 particles for two different restitution coefficient a: left a = 0.95, right a = 0.9. 
The other parameters are Ly/a — 600, Lx/a = 180 and Tt/{gLy) = 2.55 for the left panel, while Ti,/{gLy) = 5.1 for the right 
panel. 




FIG. 2: The main local observables as function of the rescaled distance y/Ly from the bottom wall: a) the local number N(y) of 
particles and the local packing fraction (t>{y) (see right y scale); b) the local mean free path X{y); c) the local rescaled granular 
temperature Tg[y)/Ti,; d) local characteristic times (see text for details). The other parameters are A*' — 1500, a = 0.95 and 
Tb/{mgLy) = 2.55. 



V^H{y)/{2v{y)) with C-ff (y) the zero-th order approximation of the homogeneous coohng rate in a layer, and 
i/(y) the kinematic viscosity in the layer; because k±{y) ex 1/X{y) this condition determines a constraint on the 
maximum number of the particles in each layer. 

• Absence of convection rolls: it is known that this granular setup is also subject to thermal convection instability 
due to the competition between gravity and temperature gradients [30, 39]; by direct inspection of the velocity 
field and the dynamical evolution of the numerical simulation we have verified that no convection is present in 
our system. 



• Fast local relaxation of microscopic modes: this condition is equivalent to ask that a particle in a layer suffers 



• 



many collisions before going out from the layer, which (in view of the fact that Li < L^) is equivalent to 
X{y) < Li (note that this is not strictly verified in the topmost layer). 

Diluteness: to avoid strong velocity correlations and hope a reasonable comparison with dilute hydrodynamic 
granular theories, we also require that the local packing fraction is small, i.e. (/>(i/) ^ 1. 

Local homogeneity: finally, a major constraint in the choice of mi is given by "quasi-homogeneity" inside the 
i-th layer, i.e. we require that [a^+i — ai]^^ai <C 1, being Ui one of the "slow" variables, i.e. N or Tg. 



After having checked by direct inspection that all these requirements are satisfied in our system and with this choice 
of TO;, we have to warn the reader about two important properties which distinguish this one with respect to other 
setups, where fluctuating hydrodynamics has been previously discussed (e.g. the homogeneous cooling state [22]). 
Each horizontal layer is considered here as a single finite system where the fluctuations of velocity shear mode is 
studied. Anyway: 

1. Each layer is an open system, exchanging particles with adjacent layers; we will discuss how this peculiarity do 
affect the studied fluctuations, in particular the autocorrelation decay. 

2. Tuning the external parameters, such as a, the total number of particles N or the parameters of the thermal 
base Th, affects the properties of each layer, e.g. local density or temperature, in a complicate and not direct 
way. Even if solutions of the hydrodynamic equations for this particular 2D setup exist [26] , the boundary layer 
near the thermal base always requires a careful treatment and makes quantitative predictions for local variables, 
from the only knowledge of external parameters, quite hard. 

IV. NUMERICAL RESULTS FOR HYDRODYNAMIC FLUCTUATIONS 
A. Time decay of the autocorrelation and dependence on wave number 

The behavior of C±{y,t) obtained in layers at different height y, from the simulations, is shown in Fig. 3a. The 
autocorrelation function of the shear mode presents two main regimes: a first rapid decay toward an intermediate 
plateau and a second decay to zero, which is roughly exponential. 

Our interpretation of the first decay is that it is associated to the openness of the layer: particles belonging to the 
layer at time t = escape from it with a typical velocity Vth{y) = \/'2Tg{y), being on average replaced by particles 
coming from the adjacent layers. In Fig. 3b we have rescaled the time with Tpxniy) — Li/vth{y), getting a fair collapse 
of the first decay in different layers. This collapse confirms that the typical time of this first decay is ~ T^xitiy) and 
that this is associated to the exchange of particles through the boundaries of the layer. The second decay can be 
reasonably fit by an exponential law. This fit works better in the central layers, i.e. 0.25 < y- < 0.65. Problems in 

the top and bottom regions could be related to the mean free path becoming too large (this is particularly true in 
the top region, where it gets closer to the width of the box): in that cases the mean free time could become close 
to the time taken by a sound wave to travel along the width of the system, producing non-exponential relaxation 
which could also be size-dependent. The region near the base is also affected by the complicate boundary layer where 
particles take an asymmetrically distributed vertical velocity because of the thermostat. From now on we focus on 
the central layers where the decay is well fit by C±_ ~ exp(— t/T(j/)). 

In those cases, it is interesting to study the decay time T{y) as a function of k and to compare it with the 
theoretical prediction associated to vorticity diffusion r^(A:, y) = l/{v{y)k'^) which describes the decay of shear modes 
autocorrelation in a homogeneous elastic system with the viscosity of the layer v{y) (see Eq. (8)). To this aim 
we have performed some simulations modifying only the horizontal size L^ of the box and maintaining the total 
density n constant. In this way we can obtain a better resolution in k. Our analysis is focalized to the total density 
ncr^ = 1.4 • 10"^ (the same of Fig. 2 and Fig. 3) and ncr^ = 1.7 • 10^^ and the comparison is shown in Fig. 4. Our first 
observation is that all data, appropriately rescaled, collapse indicating the same dependence on the local parameters 
n{y) and Tg{y) as the homogeneous system. In particular, data in the range 0.1 < fcA(y)/(27r) < 0.5 show a very good 
agreement with the prediction of simple vorticity diffusion due to shear viscosity r(y) — T^{k,y) — \/{v{y)k'^). From 
values k\{y)/{2'K) w 0.7, T{y) shows deviations from T^{y): we suspect that this is an indication of the failure of the 
hydrodynamical approach at such small lengthscales. A less clear failure is present also at very large wavelengths 
{kX(y)/(2TT) < 0.1): such a discrepancy suggests the presence of some mechanism slowing the relaxation of shear 
modes, perhaps a precursor of power-law relaxation associated to the two-dimensional geometry. 
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FIG. 3: The local rescaled correlation function C±{y, t) as function of the time t (a) and rescaled time t/Texuijj) (b) for different 
layers. The other parameters are the same of Fig. 2. 
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FIG. 4: The inverse of the decay time t{jj) as a function of the wave number k. The ordinate is been rescaled by the mean 
collision time Tc(y) = ^{y)/vth{y), while the wave number is rescaled by the local mean free path A(y)/27r. The data have a 
good collapse for the two analyzed total densities na'^ = 1.4 ■ 10~^ (open symbols) and na^ — 1.7 ■ 10~^ (closed symbols) and 
for three different layers: a) y/Ly = 0.25 (circles), b) y/Ly = 0.35 (squares) and c) y/Ly = 0.45 (diamonds). For completeness 
the theoretical prediction (dashed curve) of the local decay time Tv{y), due to the vorticity diffusion, is shown (see text for 
details). The other parameters are the same of Fig. 2. 



B. Structure factor of transverse velocity modes 

The amplitude of fluctuations of the shear mode U±{k, y, t) in a layer, which is directly proportional to the transverse 
velocity structure factor, is given (see Eq.(3)) by 



C±(fc,2/,0) = 






Niy) 



N{y}N{y) 






— ik7r{Xj —Xp)\ 



N{y) 



1 + A(fc,y) 



(17) 
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j=i p=i 
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From Eq. (17) it is seen that C±{k,y,0) can be split into two parts: a first self term which scales as ^ N{y) and a 
second term, ~ N(y)A{k, y) that contains the informations of the fc-dependence in a layer at height y. 
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FIG. 5: The local quantity A(A;, y) (see Eq. (17) for a definition) as a function of the wave number k rescaled by the local mean 
free path X{y)/{2n). The open symbols are the simulation data for the fixed layer with y/Ly = 0.35 and for two different total 
densities a) na^ = 1.4 • 10~^ (open circles) and b) na^ = 1.7 ■ 10~^ (open squares). The closed symbols are the corresponding 
values of the two densities a) and b), obtained from Eq. (18) using an exponential fit for the decay times r(j/) (data plotted in 
Fig. 4). The collapsing data are compared to the theoretical prediction (dashed curve) using in Eq. (18) the local decay time 
Tv{y) due to the vorticity diffusion. The other parameters are the same of Fig. 2. 

In a homogeneous equilibrium fluid, from equilibrium statistical mechanics, one expects the absence of spatial 
velocity correlations, corresponding to a flat C±{k,y,0), as expressed by Eq. (7). Homogeneously cooling granular 
systems (with size smaller than the critical size for linear stability of shear modes) are expected to present a non- 
flat structure factor, implying spatial velocity correlations, as marked by Eq. (12). However those systems present 
an instabihty which is not directly observed in our simulations (unless one gets much larger horizontal size L^)- 
Moreover, our layer is in a stationary state, a regime quite different from the cooling state: such a difference is 
reasonable because each layer is "driven" by particles coming from nearby layers. The best candidate as a model 
for quantitative comparison with our numerical results is the homogeneously driven system studied in [18], whose 
fluctuating hydrodynamics we have resumed in Eq. (14). Our system is different from that model, since we have not 
a uniform noise throughout the system, but the addition of energy is localized at the bottom of the box. Nevertheless 
in a single layer there is a balance between the energy exchanged with adjacent layers and the energy dissipated in 
collisions. It is therefore tempting to introduce an effective external noise whose amplitude is determined by the energy 
balance and is proportional to N {y)C,H{y)Tg{y) . Using this expression for noise and considering that the characteristic 
decay time is T(k, y), the quantity Cj_(fc, y, 0) can be written as 



C±ik,y,0) 



Niy)i, , Cniy) 



T{k,y) 



(18) 



Comparing Eqs. (17) and (18) we note that A(fc, y) can be written as Ci/(j/)T(fc, j/)/2. In an elastic homogeneous 
system this term is zero, while in a granular system we observe a non-negligible contribute of A(fc, y). In order to verify 
our assumption we have plotted in Fig. 5 the simulation data for A(fc, y) and its prediction based on the Eq. (18) and 
estimated using for T(y) the data plotted in Fig. 4. The agreement is very good for both the cases analyzed (open 
and closed symbols in Fig. 5) and this validates our considerations. Moreover we have verified that the probability 
distribution of U±{k,y,t) is close to a Gaussian, consistent with a linear Langevin equation. In Fig. 5 only results 
from the central layer {y/Ly — 0.35) are shown: however we have verified that it is equally good in the contiguous 
layers. For the sake of completeness we have shown in the figure also /S.(k,y) obtained using as characteristic time 
T^{y). The disagreement is small and analogous to that already seen in the previous section. 



C. Amplitude of fluctuations for the mode at largest wavelength 



To have a better assessment of our hypothesis (Eq. (18)), we focus now on the amplitude Cj_(fcmi„,j/,0) of the shear 
mode for the largest available wavelength kmin ■ Indeed one expects that at larger wavelengths hydrodynamics works 
better. For such a purpose, we have fixed the width L^ {Lx/(t — 180) and we have changed the total number TV of the 
particles. In this way the local density n(y) is also changed and we can explore the dependence of A(y) = A(fc,„i„, y) 
on the particle number or, equally, on 1/A(y). 
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FIG. 6: Panel a: The parameter A(y) = A(fcmin, y) rescaled by 1 
path Lx/\{jj) for different inelasticities: a — 0.95 (circles), a — 0.90 (squares) and a = 0.85 (diamonds). The dashed curve is 
the quantity C,H{y) / (2v(y)k^i„) corresponding to the homogeneous cooling theory (see Eqs. (9)-(10)) Panel b: The parameter 
Amaa; (see text for a definition), rescaled by 1 — o? , as function of the inverse of the Lx/\min, with Amin = {n"^'^^ (ij)a)~^ , 
for the same inelasticities of panel a (closed symbols). The open symbols correspond to the data obtained from Eq. (18) (see 
text for details). The dashed curve is the same of panel a valued on the layer with the maximum of the density. The other 
parameters of the two panels are the same of Fig. 2. 



Measures of lS.{y) for many different layers and for systems with three values of the restitution coefhcient a (0.85, 
0.90 and 0.95) are shown in in Fig. 6a. The data are rescaled by the factor 1 ~ a^ because this term contains the 
main dependence on a in the Eq. (18). The collapse of data is quite remarkable, if one consider that data come from 
different systems and from layers which can be at very different distance from the base: for instance, a given value of 
n{y) can correspond both to layers below the density maximum and above it. At moderate values of n{y) the data 
are closed to the values obtained using T,y(y) as decay time in Eq. (18). The disagreement at small values of n{y) is 
due to the fact that hydrodynamic hypothesis is not longer valid in this range. It is likely that the deviations from a 
clean collapse in Fig. 6a are due to the mixing of data from regions characterized by too different physical parameters 
as Tg and a. We have repeated the analysis by focusing only on the layer containing the maximum of the density 
n™"^(y), but changing the external parameters in order to explore different values of n'""^(y). The corresponding 
value Amos is presented in Fig. 6b. The data do not rescale exactly and this indicates an uncleared a dependence. 
On the other side the measures are in good agreement with the corresponding data obtained from Eq. (18) (open 
symbols in figure). This suggests that the Eq. (18 ) reproduces the right behavior provided that the suitable decay 
time is used. Concerning this point, it is important to notice that for large values of n{y)^ ^{u) (or ^max) can become 
much larger than that observed in the homogeneous cooling state (see dashed line in Fig. 6). 

V. CONCLUSIONS 

We have analyzed the fluctuations of the shear mode of the velocity in nearly homogeneous subregions of a inho- 
mogeneous granular fluid. Our conclusion is that - neglecting regions too close to the upper and lower boundaries 
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- the large time temporal decay of their autocorrelation does not deviate dramatically from the expectation for a 
homogeneous equilibrium fluid, i.e. \/v{y)k'^ , provided that hydrodynamic (large enough) scales are considered. At 
too large scales, comparable with the total horizontal size, some deviations from the \/v{y)k'^ behavior is appreciated, 
likely due to finite size effects and problems associated to the 2D geometry. 

The study of the amplitude of fluctuations, i.e. the transverse velocity structure factor, suggests that for the central 
layers, the effective noise is well described by a sum of two contributions, both white and Gaussian: an internal one, 
associated to the vorticity current, which is at "equilibrium" with the local temperature, and an external one which 
is responsible for the balance of local energy which is continuously lost in collisions. We have shown that those as- 
sumptions lead to a good estimate of the structure factor in all cases analyzed, if the autocorrelation decay time used 
in the formula is the one measured in the simulations: such an agreement holds even when the decay time deviates 
from the simple prediction l/v(y)k'^ . 

We stress that formula (18) for the amplitude of fluctuations represents an explicit violation of the equilibrium 
Fluctuation-Dissipation Relation, which would yield a constant structure factor: such a violation can be quite large 
with respect to those discussed in the homogeneous cooling regime [22]. In particular, violations observed here lead 
to an amplitude of fluctuations which can be much larger than that expected at equilibrium or in the homogeneous 
cooling theory. It is also remarkable to notice that the correlations observed here are not related to the gradients in 
the system (as it happens for other correlations in molecular systems with an imposed thermal or density gradient), 
but are a direct consequence of inelastic collisions. 
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